Vortex nucleation and hysteresis phenomena in rotating Bose-Einstein condensates. 
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We study the generation of vortices in rotating Bose-Einstein condensates, a situation which has 
been realized in a recent experiment (K. W. Madison, F. Chevy, W. Wohlleben, J. Dalibard, Phys. 
Rev. Lett. 84 806 (2000)). By combining a linear stability analysis with the global optimization 
of the nonlinear free energy functional, we study the regimes that can be reached in current exper- 
iments. We find a hysteresis phenomenon in the vortex nucleation due to the metastabilization of 
the vortexless condensate. We also prove that for a fast enough rotating trap the ground state of 
the condensate hosts one or more bent vortex lines. 

PACS number(s): 03.75.Fi, 05.30. Jp, 67.57.De, 67.57.Fg 



Vortices, or vortex lines we should rather say, con- 
stitute the most relevant topological defect in Physics. 
They consist on a twist of the phase of a wave function 
around an open line and they are typically associated to 
a rotation of a fluid whatever the fluid is made of (real 
fluids, optical fluids, quantum fluids, ...) Q. Vortices 
are one of the means by which quantum systems acquire 
angular momentum and react to perturbations of the en- 
vironment. They have already been predicted, observed 
and studied in the superfluid phase of 4 He and are in- 
deed known to be the key to some important processes 
in these systems, such as dissipation, moments of inertia 
and breakdown of superfluidity. This is why extensive re- 
search on vortex generation, stability and dynamics has 
been conducted in the field of Bose-Einstein condensation 
(BEC) in the last years 

Vortices and other defects usually involve more energy 
than other equilibrium states, such as convex nodeless 
states. Therefore, in order to produce a vortex in a 
condensate, one must induce a transition from a uncon- 
densed or condensed convex cloud to the desired state, 
by means of an external action such as a change of the 
confining potential. 

For instance, in recent experiments performed by a 
group from the Ecole Normale Superieure (ENS) ||, vor- 
tices are created by rotating an elongated trap which is 
slightlty deformed along its transverse dimensions. After 
this preparation, the trap is switched off, and the conden- 
sate expands until vortices become directly observable. 

In these experiments there are several controversial 
points which are: (i) Vortices first nucleate at a trap 
rotation speed or critical frequency, Q\, larger than that 
required to stabilize a vortex line in the Thomas-Fermi 
theory, (ii) when seen from the top vortex cores seem 
partially filled and (iii) after the nucleation of the first 
vortex the angular momentum grows continuously with 
the rotation speed (not only with discontinuous jumps). 

In this paper we reach a global view of the transitions 
between equilibrium states which are induced by the trap 
rotation, plus a simple explanation of the most contro- 



versial points. Our study consists of two parts. First we 
focus on the analysis of stationary states and discuss the 
stability properties of the simplest solutions, i.e. straight 
vortices, which allows us to get some insight on the prob- 
lem. Our main result is that the metastabilization of a 
vortexless state induces hysteresis in the vortex nucle- 
ation process. In other words, the trap rotation speed 
must exceed a critical value, CIm > to produce vortices, but 
that speed is well beyond the values, fii,f22, ■ ■ •, which 
are required to make vortices energetically favorable. In 
the second part we find numerically the actual ground 
states of the condensate for different angular speeds. We 
will show that once the first vortex is nucleated the angu- 
lar momentum grows almost continuously with the speed 
of the trap, while the ground state (GS) mutates into dif- 
ferent deformed states. 

Several other works have contributed to the theoreti- 
cal description of the experimental "anomalies". First, 
Ref. in Q the dynamics of vortex lines is studied ana- 
lytically, and different modes are obtained which reflect 
the transverse tension of vortex lines. Next, in Ref. || 
the existence of these modes is confirmed numerically 
a relation between these modes and the large value of 
the critical angular speed of the ENS experiments is pro- 
posed. Finally in Ref. || the authors derive a condition 
for the efficiency of external perturbations when trying 
to induce a mechanical response in the condensate. 

The model.- For most of current experiments it is an 
accurate approximation to use a zero temperature many- 
body theory of the condensate. In that limit the whole 
condensate is described by a single wavefunction ip(r,t) 
ruled by a Gross-Pitaevskii equation (GPE) [ [Tcfl . 

In Ref. j^] the trap is initially harmonic with axial sym- 
metry, but then a laser is applied which deforms it and 
makes it rotate with uniform angular speed £1. On the 
mobile reference frame which rotates with the trap, the 
experiment is modeled by a modified GPE 
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Here L z = i (xd y — yd x ) is the hermitian operator that 
represents the angular momentum along the z-axis and 
the effective trapping potential in given by 

Voir) = - e). 2 + + s)/ + l r Wz\ (2) 

In Eq. (Q) we have applied a convenient adimen- 
sionalization which uses the harmonic oscillator length, 
a± = y/h/mRi ) u)± and period, r = wj 1 . With these units 
the nonlinear parameter becomes g = 47ras/aj_, where 
as — 5.5 nm is the scattering length for 87 Rb, the gas 
used in Ref. p). Following the experiment we will take 
oj z = 2tc x 11.6 Hz and ui± = 2ir x 232 Hz, and we will 
use Ng = 10000, which corresponds to a few times 10 5 
Rb atoms. For the small transverse deformation of the 
trap we have tried e = 0.0, 0.03 and 0.06 (e = 0.03 is the 
closest one to the actual experiment ||). 

The norm, N[tp] = J M dr, which is related to the 
number of bosons in the condensate, and the energy 



We must note that the condition £1 > Jl m , does 
not imply that ip m is a global minimum but only that 
E(ipo) > E(ip m ). It is feasible for other states with the 
same vorticity to have less energy than the straight vor- 
tices. That is, another state ij) m might exist such that 
E(4>m) < E(ip m ). This is known to be the case for ip2 
which is energetically less favorable than a pair of m = 1 
vortices §,[§. 

In fact, for an axisymmetric vortex ijj m to be stable, it 
should be at least a local minimum of E(ip). One way 
of studying the local stability is by using linear stability 
analysis, i. e. we linearize the energy functional around 
any state <5(r) 

Ety + S] = E[iJj] + E'[^}{5) + E"[iJ;}(6,6) + 0(\5\ 3 ) (4) 



and find a quadratic expansion 



E"[i>}(6,5) = J2*m,Mi>M 
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are conserved quantities. Each stationary solution of the 
GPE of the form ipjr, t) 
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in which the perturbation is expanded on a suitable ba- 
sis mj (p,z)e™ e | as S(r) = E m Jm^ m .j(P^y me - 
Following this procedure we are able to study the curva- 
tures of the surface (tp, E[tjj]) around that state — namely, 
r), is a critical point of tne eigenvalues X m>j (Q;ip) = A mj (0;V) - mfi— . This 
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the energy (|[) with the norm constraint, 

Critical frequencies.- Each state of the system is a 
"point" in the infinite-dimensional functional space on 
which the energy Eo(ip) and the angular momentum 
L z (ip) are defined. In principle, due to the conserva- 
tion of the energy, and due to the fact that dissipation 
at most drains energy out of the system, the minima of 
Eo(ip) are associated to stable states. Around these min- 
ima, the "motion" of the system is confined by energy 
barriers. However, rotation involves a pointwise change 
of the height of these infinite dimensional "surfaces" of 
energy, which move from Eo(tp) to Eo(tp) — QL z (ip). This 
shift, which is stationary for each dynamical configura- 
tion of the trap adopted, can turn stable local minima 
into saddle points and open paths for the evolution of 
the condensate from simple nodeless state to states with 
one or more vortices. 

Let us consider the case of an axially symmetric trap. 
These traps admit axially symmetric solutions with ei- 
ther no vortices, ?po, or with one centered vortex of 
integer charge, ?/> m oc expjim arctan(y/a;)}. In non- 
rotating traps, the vortexless state is actually the one 
with less energy. As the rotation speed is increased, 
the energies of the states with positive vorticity are all 
shifted down, E rn (Q) — E rn (0) — mf!, so that when 
fl > Q rn = E ™(°)^ E °( n ) ; m = l } 2, the state ipm, i- e. 
a symmetric and straight vortex with vorticity m, is en- 
ergetically more favorable than tpQ. The values Vl m are 
thus first estimates for the critical frequencies at which 
states with vorticity m could be created. 



analysis allows us to decide when the operator E" [ip m ] is 
positive definite and consequently when a particular ip m 
becomes a local minimum. This happens for a critical 
frequency that will be denoted by 
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For instance, when £1 — Cli the straight vortex becomes 
a local minimum. 

To generate a condensate with vorticity two conditions 
must be satisfied. It is clear that the trap must rotate 
fast enough to make a vortex state energetically favor- 
able, but it is also necessary that an energy decreasing 
path exists from ?/>o to the vortex state, otherwise the 
transition could be inhibited due to the energy barrier 
around i/'o making this state metastable. Mathemati- 
cally there should exist a perturbation (a "direction" in 
the phase space), 4> m ^, such that £1 > \ mt i(0;i/jo)/m, a 
condition which looks similar to the bounds of Ref. [||. 
We will define a destabilization frequency, given by 



Qm = rnin 

m,j 
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Within the interval fl e [0, Q.m) a condensate in a vor- 
texless state cannot be forced to acquire angular momen- 
tum unless some additional energy is pumped in and the 
nodeless state i/jq could be metastable. 

Results for symmetric vortices.- One could be tempted 
to think that the real situation is the simplest one, i. e. 
that by increasing Q one should reach a point Qi in which 



2 



E(ipo) > E(ipi) and at the same time fii > CIi,Hm- In 
this situation ip\ would become also a local minimum 
(Oi) and the core state would loose its stability (Om)- 
However this is not the case, and the situation is more 
complicated as we will see below. 

With the previous definitions in mind we have com- 
puted numerically the states tpQ, ipi, ... as well as the dif- 
ferent frequencies Q%, fi2, O3, ^1 and Om- The results are 
plotted in Fig. [I] for a very elongated trap such as those 
used in the ENS experiment 0. 

There are two relevant conclusions that may be ob- 
tained from this picture. The first one is that Ox > 
Oi , which means that the ground state of the sys- 
tem may never be a symmetric vortex line. In other 
words, should it be energetically favorable for the ground 
state to acquire some vorticity, it will never be by means 
of straight vortex lines, but with some other structure. 
From a practical point of view, the large difference be- 
tween fii and fii also implies that the ENS experiment 
is working on a regime in which straight vortex lines are 
very hard to obtain. What are the new vortex structures 
which appear instead above £l\ is a point we will discuss 
later. 
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FIG. 1. (a) Critical angular speeds for an elon- 

gated condensate (uj z /lu± — 18.7) at different nonlinearities 
(U — gN = 4irNa Rb /aj_). We plot the Thomas- Fermi esti- 
mate for the energy difference between a core and vortex state 
(Qtf, dash-dot line) the speeds fh, ^2 and at which vor- 
tices with charge m — 1,2,3 become energetically favorable 
(solid lines), and the angular speed at which ipo (vortexless 
state) becomes a saddle point (!1m, dashed line), (b) Regions 
with different phenomenology: from A to C there is an energy 
barrier surrounding the vortexless state; in A, B and C the 
ground state has vorticity equal to 0, 1 and 2 respectively; in 
D the energy barrier around the core state vanishes and more 
and more vortices become feasible. 

The second and most relevant feature is that the vor- 
texless state, ipoi remains a local minimum up to a speed, 
f^M) which is much larger than the value at which vor- 
tices become energetically favorable (Regions B and C in 
Fig. 0(b)). Before reaching £1 = Qjm it is energetically 
expensive to introduce angular momentum into the con- 
densate, since an energy barrier must be surpassed. For 



fl > f2j\/ states with one or more vortices become fea- 
sible at the same time. This persistent metastability of 
tpo should be responsible for the experimental observa- 
tion that a single angular speed is capable of producing 
states with a different number of vortex lines Q . 

Our findings imply that once a vortex state is reached, 
the rotation speed may be ramped down and vortices 
should remain stable for some range of Q values. For 
instance, if after reaching Qjf ~ 0.63uj± for U = gN = 
3000 one gets a state with a vortex, then the rotation 
speed may be ramped down until fl\ ~ 0.46wj_ giving rise 
to a hysteresis phenomenon as it is graphically illustrated 
in Fig. |(b). 




FIG. 2. (a) Angular momentum, L z , of the ground state 
of the energy functional E(i/j) [Eq. (g)] as a function of Q. 
Circles correspond to solutions over a grid with 32 x 32 x 64 
Fourier modes for s = 0, 0.03, 0.06. The solid line corresponds 
to the solutions on a grid with 64 x 64 x 128 modes and 
e = 0.03. (b) Schematic picture of the angular momentum 
of the ground state versus angular speed, with a graphical 
description of the hysteresis mechanism. All figures are adi- 
mensional. 

True ground states.- For f2 > Qmi or if the the energy 
barrier which surrounds the core state is overstepped, 
a condensate produced by means of evaporative cooling 
should correspond to the absolute minimum of the en- 
ergy in the configuration space. We have worked nu- 
merically with the energy functional [Eq. in three 
spatial dimensions using a technique known as Sobolev 
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gradients to find the ground states subject to some rea- 
sonable constraints — i. e. the norm and angular speed — 
(simple minimization methods do not work for this prob- 
lem) . The details of the procedure are given in Ref. . 
We have applied this method on a Fourier basis with 
32 x 32 x 64 modes (enough for plotting purposes) and 
with 64 x 64 x 128 modes (which is needed to lower the 
error in L z below 1%.) 

Figs. ^ and || summarize our results. In Fig. ||(a) 
we see that the ground state acquires vorticity when Q, 
is smaller than the core-state destabilization speed and 
the experimental values found in Q, in agreement with 
the predictions of the preceding paragraphs. Our numer- 
ical method also solves the question posed before: if the 
ground state must have some vorticity and it is not a 
straight vortex, what is its shape? In Fig. || we see that 
vortices are nucleated with deformed shapes. The com- 
bination of the sudden destabilization of the vortexless 
state, and the longitudinal deformations of vortex lines 
is enough to explain why vortex lines seem partially filled 
and why the angular momentum evolves almost contin- 
uously above f2j\/ ■ More details on the structure of bent 
vortices will be given elsewhere 111] . 




FIG. 3. Shape of the ground state of the Hamiltonian 
(|) for Q = 0.75cj ± ,7 = 18.7, Ng = 1000 and e = 0.06. 
(a) Transverse section, (b) Density plot. One quarter of the 
condensate has been removed to allow the direct observation 
of the inner structure of the vortex line. 

Discussion and conclusions.- We have used the GPE 
to study the nucleation of vortices in an elongated trap 
similar to that used at the ENS || . Our main prediction 
is that in such elongated traps there exists a mechanism 
which prevents the nucleation of vortices unless the con- 
densate rotates much faster than the speed required to 
make one or more vortices energetically stable (metasta- 
bilization of ipo). 

We predict that in such elongated traps, straight vor- 
tex lines are locally stable only for very large angular 
speeds — well beyond the values at which higher vor- 
ticites become preferable — . And we find the longitudinal 
deformation of these topological defects to be responsible 
for an almost continuous growth of the angular momen- 



tum with respect to the angular speed. 

It is difficult to observe the bending of vortex lines 
in current experimental setups since condensates are ex- 
panded too much along their transverse dimensions. Nev- 
ertheless, other predictions of this paper such as the 
metastabilization phenomenon should be easily testable. 
First, before the destabilization of the vortexless state, 
f2jU) the vortex nucleation process is prevented, and the 
condensate may only adapt to the rotation speed by 
means of transverse deformations which are described in 
p3| . And second, a more controlled nucleation of vor- 
tices is possible by ramping up the trap beyond SIm and 
then carefully slowing down the condensate to some point 
above Qi . An expansion of such cloud should lead to the 
observation of vortices for much lower values of fl than 
those reported in Ref. ||. 
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[1] P. G. Saffman, "Vortex dynamics" , Cambridge University 
Press (1997); Y. Kivshar, B. Luther-Davis, Phys. Rep. 
298, 81 (1998). 

[2] D. S. Rokhsar, Phys. Rev. Lett. 79 (1997) 2164; R. J. 
Dodd et al, Phys. Rev. A 56 (1997) 587; R. Dum, et 
ai, Phys. Rev. Lett. 80 (1998) 2972; F. Zambelli, S. 
Stringari, Phys. Rev. Lett. 81 (1998) 1754; M. Caradoc- 
Davies et ai, Phys. Rev. Lett 83, 895 (1999) T. Isoshima 
and K. Machida, Phys. Rev. A 59 (1999) 2203; D. L. 
Feder et al, Phys. Rev. Lett. 82 4956 (1999); M. R. 
Matthews et al, Phys. Rev. Lett. 83 2498 (1999); D. L. 
Feder et ai, Phys. Rev. A 61 011601 (2000); A. Svidzin- 
sky and A. L. Fetter, Phys. Rev. Lett. 84, 5919 (2000). 

[3] D. A. Butts and D. S. Rokhsar, Nature 397, 327 (1999). 

[4] J. J. Garcfa-Ripoll and V. M. Perez-Garcfa, Phys. Rev. 
A 60, 4864 (1999). 

[5] M. R. Matthews, B. P. Anderson, P. C. Haljan, C. E. 
Wiemann, E. A. Cornell, Phys. Rev. Lett. 83 (1999) 2498; 
J. J. Garcfa-Ripoll, V. M. Perez-Garcfa, Phys. Rev. Lett. 
84, 4264 (2000); V. M. Perez-Garcfa, J. J. Garcfa-Ripoll, 
Phys. Rev. A 62 033601 (2000). 

[6] K. W. Madison, F. Chevy, W. Wohfleben, J. Dalibard, 
Phys. Rev. Lett. 84 806 (2000); F. Chevy, K. W. Madi- 
son, J. Dalibard, Phys. Rev. Lett. 85, 2223 (2000). 

[7] A. A. Svidzinsky, A. L. Fetter, Phys. Rev. A, 62, 063617 
(2000). 

[8] D. L. Feder, A. A Svidzinsky, A. L. Fetter, C. W. Clark, 

|cond-mat/0009086| 
[9] F. Dalfovo, S. Stringari, e-print |cond-mat/0009187 
[10] F. Dalfovo et ai, Rev. Mod. Phys. 71, 463 (1999). 
[11] J. J. Garcfa-Ripoll, V. M. Perez-G arcfa, e-print 

[http://xxx.lanl.gov/abs/mat.AP/0008225l . 
[12] J. J. Garcfa-Ripoll and V. M. Perez-Garcfa (preprint). 
[13] J. J. Garcfa-Ripoll and V. M. Perez-Garcfa, e-print cond- 
mat/003451. 



4 



